Stochastic dynamics of a Josephson junction threshold detector 
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We generalize the stochastic path integral formalism by considering Hamiltonian dynamics in the 
presence of general Markovian noise. Kramers' solution of the activation rate for escape over a barrier 
is generalized for non-Gaussian driving noise in both the overdamped and underdamped limit. We 
apply our general results to a Josephson junction detector measuring the electron counting statistics 
of a mesoscopic conductor. Activation rate dependence on the third current cumulant includes an 
additional term originating from the back-action of the measurement circuit. 
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I. INTRODUCTION 

Detecting electron counting statistics has become a 
major experimental challenge in mesoscopic physics. 
First attempts to measure non-Gaussian effects in current 
noise have revealed that the detection problem is quite 
subtle. In particular, the experiment 1 found that the 
third current cumulant was not described by the simple 
theoretical prediction, 2 but was masked by the influence 
of the measurement circuit causing an additional "cas- 
cade" correction. 3 ' 4 Recent experiments demonstrated a 
measurement of the third current cumulant without cas- 
cade corrections, and the detection of individual electron 
counting statistics. 6 Stringent bandwidth requirements 
in measuring the third cumulant suggested that further 
experimental advances would require a new approach. 

A conceptually different way to measure rare current 
fluctuations is with a threshold detector, 7,8 the basic 
idea of which is analogous to a pole vault: A detection 
event occurs when the measured system variable exceeds 
a given value. A natural candidate for such a detector is a 
metastable system operating on an activation principle. 9 
By measuring the rate of switching out of the metastable 
state, information about the statistical properties of the 
noise driving the system may be extracted. A threshold 
detector using an on-chip conductor which contains a re- 
gion of negative differential resistance 8 ' 10 was proposed 
by the authors and shown to be capable of measuring 
large deviations of current. Tobiska and Nazarov pro- 
posed a Josephson junction (J J) threshold detector, 7 the 
simplest variant of which (see Fig. 1) operates essentially 
in a Gaussian regime. 11 The third cumulant contribution 
is small 12 and may be extracted using the asymmetry of 
the switching rate with respect to bias current, as has 
been demonstrated in recent experiments. 13,14 

In this Letter we solve Kramers' problem 9 of noise- 
activated escape from a metastable state beyond the 
Gaussian noise approximation and investigate how the 
measurement circuit affects threshold detection. Start- 
ing with general Hamiltonian-Langevin equations which 
includes deterministic dynamics, dissipation, and fluctu- 
ations, we represent the solution as a stochastic path in- 
tegral of Hamiltonian form 15,16 by doubling the number 
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FIG. 1: Josephson junction (JJ) threshold detector: Simpli- 
fied electrical circuit with the JJ (marked with an X) and 
mesoscopic system inserted. The driving noise of the system 
creates fluctuations of the center node voltage V, that can 
activate the JJ out of its supercurrent state, into its running 
dissipative state until it is recaptured. 

of degrees of freedom. In the weak damping case, the 
dynamics is dominated by energy diffusion, which we ac- 
count for by a change of variables, enabling an effectively 
two-dimensional representation. We calculate the escape 
rate via an instanton calculation, and obtain a formal 
solution of Kramers' problem. 9 Applying these general 
results to a JJ threshold detector, we account for the in- 
fluence of the measurement circuit and find that the cas- 
cade corrections are a consequence of the non-equilibrium 
character of the noise. 17 



II. HAMILTONIAN-LANGEVIN EQUATIONS 

In order to reformulate and solve Kramers' problem 
beyond Gaussian noise, we introduce a theoretical frame- 
work that relies on a separation of time scales: Slow 
motion of a deterministic system on a time scale Tq is 
affected by quickly fluctuating noise sources with corre- 
lation time tq <C Tq. Quite generally, this situation can 
be described by Hamiltonian-Langevin (HL) equations. 
For instance, in a two-dimensional phase space (p, q) the 
equations of motion are 

q = dH(p,q)/dp + I g , p=-dH(p,q)/dq + I p , (1) 

where H(p, q) is the Hamiltonian that generates deter- 
ministic motion, and I q and I p are white noise sources. 



2 



On the intermediate time scale t, such that To 3> t 3> To, 
the sources are Markovian and fully characterized by the 
generating function H.(X p ,X q ) of their cumulants (irre- 
ducible moments): ((1^1™)) = d^ p d^H(0,0). 

While Langevin equations may be solved by stan- 
dard methods in stochastic physics, 18 the methods fail 
when applied to higher-order cumulants. The reason is 
that slow variations of p and q affect the noise sources 
and lead to cascade corrections, 3 ' 16 which cannot be ob- 
tained from linearized differential equations (1). Follow- 
ing the steps of Rcfs. [15,16], we represent the slow evo- 
lution of the system by the stochastic path integral (SPI) 
P = JT>A jDRexp(5), where the action S is given in the 
explicitly canonically invariant form as 

S = J dt'[—A • R + A • {R, H} + H(A, R)]. (2) 

Here R = (p, q) and A = (X p , X q ) are the sets of physical 
and canonically conjugated "counting" variables, respec- 
tively, and {• • •} denotes the Poisson bracket with respect 
to p and q. By fixing R in the final state of (2) we obtain 
the probability distribution T(R), while fixing the final 
A variables turns the SPI into the moment generating 
function -P(A). 

The large parameter T /t S> 1 allows the saddle- 
point evaluation of the SPI and thus requires solving 
Hamilton's equations of motion in the extended space. 16 
There always exists a trivial solution A = and R = 
{R, H} + (I), where I = (I p , I q ), that describes the "aver- 
age" dynamics in physical space with a null action S = 0, 
giving the proper normalization of the distribution P. In 
the context of the generalized Kramers' problem, one has 
to find a non- trivial instanton solution A; n (R) of Eq. (2), 
for which the SPI gives a rate of the noise activated es- 
cape from a localized state, T cx exp (£;„)■ 19 

III. ENERGY DIFFUSION AND ESCAPE RATE 

As an important example, let us consider quasi- 
periodic motion with the period To 3> tq. We change 
variables to energy E — H(p, q) and "time" s = s(p, q), 20 
accompanied by a new set of counting variables (X B , A s ), 
which are defined via X p = Xe9 p H + X s d p s and X q = 
\ E d q H + \ s d q s. Using {E,H} = and {s,H} = 1, we 
obtain A • {R, H} = X s . Fluctuations of the s variable 
increase the action without leading to escape, therefore 
they can be neglected by choosing A s = 0. In the weak 
damping regime, one can set s = t, so that X p = C[\e and 
\ q = —pX B - The action for the energy diffusion then 
reads, 

S = J dt'[-X E E + H(X E q,-X E p)}. (3) 

The following intuitive argument supports this result. 
The energy balance equation E = ql p — pl q follows di- 
rectly from the HL equations (1) and takes the form of 



a Langevin equation. If To 3> tq, then variables q and p 
change slowly and can be considered as being "effective 
charges", which after following [16] leads to Eq. (3). 

To leading order in weak damping we can replace the 
generator H in Eq. (3) with its average over the oscil- 
lation period (H)e = {^/Tq) § dtli, evaluated for fixed 
X E and E. Corrections in damping may be found by 
taking into account slow energy dissipation E = d\ E H 
and Xe = —d E T~t, while averaging over the period T . 
We are interested in the instanton solution Xe = Ai n (-E) 
with (H)e = in the initial and final state. 10 Since the 
"Hamiltonian" (Tt) e is an integral of motion, we obtain 
an important result for the escape rate, 

logT = - J X in dE, (H(X in q, -X inP )) E = 0, (4) 

which formally solves Kramers' problem for arbitrary 
Markovian noise in the weak damping limit. Below we 
apply the theory outlined here to the stochastic dynamics 
of a Josephson threshold detector. 

IV. JOSEPHSON THRESHOLD DETECTOR 

The circuit in Fig. 1 shows the essential part of the 
detector comprised of the J J with Josephson energy Ej, 
and the capacitor, C. The circuit is current-biased with 
Ib through the macroscopic conductor and by the system 
current Is, which is to be measured. According to Kirch- 
hoff 's law, the total bias current Is + Ib is equal to the 
sum of the Josephson current 1$ = (_E//<f> ) sin0 where 
<I>o = h/2e, and the displacement current Iq = CV. This 
leads to the equation of motion for the superconducting 
phase (f>, 

C<S> 2 o 4> + Ejsm0 = <t> o (Is + I B ), (5) 

where we used the relation V — $o0- 

In order to simplify the following analysis and concen- 
trate on our main message, we make a number of as- 
sumptions, most of which will be relaxed later. First, 
we consider an ohmic system and bias resistor, so that 
(Is) = Js~ G S V and (I B ) = J B - G B V, where G s , G B 
are the system and bias conductances, and the constant 
currents Js — GsVs, J B = G B V B are just tunable pa- 
rameters. The bias resistor, being a macroscopic system, 
creates Gaussian Nyquist noise {(I B )) = 2TG B . We fur- 
ther assume a high-impedance circuit, so that the back 
flow part of the bias current G B V and the Nyquist noise 
may be neglected, I B = J B , which we refer to as the ideal 
detection scheme. Then Eq. (5) can be rewritten as a set 
of HL equations (1) for the phase variable (j> and canon- 
ically conjugated momentum p = (fr Q (where Q = CV 
is the total charge on the capacitor) 

0=p/m, p=-dU/d<f> + ^ (I s -Js), (6) 
with "mass" m = ^C. 
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Equations (6) describe the motion of a "particle" in 
the tilted periodic potential 



U{4>) = -Ej cos (j> - $ ( J s + Je 



(7) 



stimulated by the dissipative part Ig — Jg of the system's 
current. For later convenience, we define the normalized 
total bias current as J — &o(Js + Jb)/Ej. Dissipation 
leads to relaxation of the JJ into one of its metastable 
states. In the supercurrent state, the phase is localized 
in one of the potential wells, so that (V) = 0. In the 
dissipative state, the phase drifts along the bias which 
generates the non-zero voltage drop V. In the remaining 
part of the paper we investigate noise activated escape 
from the localized state. 



V. WEAK DAMPING REGIME 

Here the system conductance is small, Gg <C ui p \C, 
so the phase oscillates with the plasma frequency u> p i — 
- J) 1/4 , where ilj = (Ej/^C) 1 / 2 . The energy 
relaxes slowly with the rate Gg/C to the local poten- 
tial minimum. Wc further assume the separation of time 
scales, 1/tq ~ max{eVs,T} 3> ftiv p i, so that the noise 
source Ig is Markovian. Comparing (6) and (1), we iden- 



tify q 



I P = &o{I,s — Js) and I q = 0, so the equations 



for the escape rate and "instanton line" read 



io g r 



X in dE, (H(\ ia <P)) B = 0, (8) 



where TL is the generator of the cumulants of the dissi- 
pative part of the system current, Ig — Jg, and <&o plays 
the role of an effective charge. 

We first consider the system in thermal equilibrium 
with Gaussian Nyquist noise, ((Ig)) = 2TGg, as a check- 
point of the theory. The cumulant generator acquires the 
simple form H = -$ Gg4>(\ in <j)) + TGg(\ in (j)) 2 , where 
the first term comes from the linear response of the sys- 
tem current Ig — Jg = —GgV to the potential V = &o<fi, 
while the second term is the noise contribution. Averag- 
ing 7i over the period of oscillations, we observe that the 
term ((f) 2 ) e cancels, so that A;„ = $o/T. Using Eq. (8) 
we obtain Kramers' well-known formula for the rate of 
thermally activated escape 



io g r = -Ac//r, 



(9) 



where the potential threshold AU is a function of nor- 
malized bias J. For the potential (7) one obtains 
AU/Ej = 2(1 - J 2 ) 1 / 2 - 2JarccosJ r , see Fig. 2. Lim- 
iting values are AU/Ej = 2 for J — 0, while as J — > 1, 
AU/Ej^(4V2/3)(l-J) 3 / 2 . 

The above example shows that the argument of 7i 
is small. Indeed, we estimate <j> ~ u> p \, so e</>Aj n ~ 
fku p i/T <C 1 due to the separation of time scales. There- 
fore, for any system away from equilibrium we can ex- 
pand (7€)e in Eq. (8), average the result over the period 
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FIG. 2: The normalized potential threshold AU/Ej and di- 
mensionless factors D\ and D2, are plotted versus the nor- 
malized current bias J = $o(^s + Jb)/Ej. Inset shows the 
phase space plot of a trajectory for weak damping (thick line) 
which leads to the escape, and the separatrix (dashed line). 
The turning points (j>o are shown by the dots. 

of oscillations, and invert the series for A; n . To second or- 
der in Ai n we again obtain Kramers' formula (9) with the 
temperature T replaced with the effective temperature 



T, 



eff 



«4»/2Gc 



(10) 



of the non-equilibrium noise. 

To find the third cumulant correction to the Kramers' 
formula, we expand (W(0A; n )) e in Ai„ and the voltage 
V, and collect all terms to next order in small parameter 
fvup\/T c ff. Two terms which contain ((Ig)) and d\/ s ((Ig)) 
come with the factor (4> 3 )e 21 This factor vanishes to 
leading order in damping, because for fixed energy <j) 
is an odd function of time over the oscillation period. 
Therefore, we evaluate the average (<P 3 )e by taking into 
account slow variation of the energy. This is accom- 
plished by taking functional variations of the action to 
find A in = $ /T e ff and E = § 2 Gg(j) 2 . We skip a number 
of straightforward steps (see Appendix) and present the 
result for the non-equilibrium escape rate 



io g r 



AU -Di$o£./((4»tot 



dE ((0o 
Ej 3 



)E 



(11a) 
(lib) 



where the dimensionless factor D\(J) is a characteristic 
of the detector. Appearing in the integral (lib) is one 
of the turning points of the oscillating phase, <po, defined 
by the solution of U(4>q) = E that is nearest the top of 
the potential. For the potential (7) -Di(O) = 7r/3 and for 
strong bias J -> 1 we obtain D X (J) w 0.8(1 - J) 2 . The 
parameter D\, evaluated numerically, is shown in Fig. 2. 
The total third current cumulant 



((4hot 



am 



3T eS d Vs ((I 2 )) 



(12) 
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is taken at V = and contains a correction originating 
from the slowly varying second current cumulant. This 
contribution is analogous to the cascade correction 3,4 di- 
rectly observed in third cumulant. 1 The correction is gen- 
erally not small and may even change the sign of the total 
cumulant. For instance, for systems far from equilibrium 
((1^)) — F n (Is), where F n are the dimensionless Fano 
factors, we obtain ((i|)) to t = (F 3 - 3F 2 /2)(I S ). 



VI. OVERDAMPED REGIME 

In this regime the conductance is large, G5 » w pl G, 
and the dynamics is entirely due to slow phase relaxation 
with the rate ufyC/Gs- Therefore, we can set p = m<fi 
and A<j = 0, and neglect the first term in the action (2), 
so that the action reads: S = J dt [W($ A p , 4>) — X p d^U], 
where again H is the generator of the cumulants of the 
dissipative part of the system current, Is — Js- 

The following analysis is analogous to that of the 
weak damping regime. We first expand TL to second or- 
der in X p and find from the equations of motion that 

4> = T cS X p and X p = d^U / (<5>lT cS G s )- Substituting 
these results back to the action we immediately obtain 
Kramers' formula (9) with T replaced with T e ff. Next 
we observe that the argument of H is small, namely 
e&oXp ~ huipiC/GsTeff <C 1 due to the separation of 
time scales requirement. We collect all the terms to third 
order in this parameter and evaluate them perturbatively 
using the above results for A p and <j). We finally obtain 



logT = - — 



AU D 2 Ej((ll)) tot 
T cS 4>oT 3 ff G 2 ' 



Do = 



(1/6) J < 



>(d^U/Ejf 



(13a) 
(13b) 



where the total cumulant ((/|))tot is given by Eq. (12), 
and D2 (J) is a dimensionless detector property. For the 
potential (7) it is given by D 2 = (1/6) {1+2 J 2 ) arccos,/- 
(1/2)J(1 - J 2 ) 1/2 , see Fig. 2. Limits are D 2 (0) = tt/12, 
and as J -» 1, D 2 {J) w 0.25(1 - Jf' 2 . 



VII. DISCUSSION 

We now remark on the application of our results to 
the detection of non-Gaussian fluctuations. It is evident 
from Fig. 2 that D lfi < MJ/Ej as J -> 1. Therefore, in 
the strong bias regime the third cumulant contribution 
in (11) and (13) is suppressed compared to the Kramers' 
term So = AU/T c s. But even for a relatively weak bias, 
when D\ t 2 ~ 1, non-Gaussian effects are small. Indeed, it 
is easy to estimate the correction as (7£ w d/Q)£o, where 
the JJ quality factor Q = Cco p i/Gs > 1 in the weak 
damping regime (wd), and the ratio lZ w d = hjJ v \/T e R < 1 
due to the separation of time scales. Similarly, in the 
strong damping regime (sd), Q < 1, the correction is of 
order 1Z so \So, where lZ s d = huj 2 y C '/GsT c g < 1 due to the 



separation of time scales. Since Sq itself cannot be too 
large for the escape to be detected, in experiments one 
should try to saturate the above inequalities and use the 
asymmetry of the third cumulant (12) as a function of 
the current bias Jb- 

Next, we briefly discuss nonideal detection and non- 
linear effects. The finite conductance of the bias resistor 
Gb contributes to the total conductance of the circuit 
Gtot = G s + G B and the Nyquist noise ((1%)) = 2TG B 
adds to the total noise power. Therefore, in Eqs. (11- 
13) one has to replace Gs with G to t and the effective 
temperature with T = (GgT e fj + GsT)/Gtot- The non- 
linearity of the system current leads to the correction 
(X p /2)(^ (j)) 2 dy s (Is) to H, and thereby to an additional 
contribution to the total third cumulant, ((/|))tot = 
((4)) - 3fd Vs ((Is)) + 3f ,2 ^ s (7 s ), in both transport 
regimes considered above. Quite remarkably, this total 
third cumulant is related via ((/|))tot = 3G 2 Gtot((^ 3 )) to 
the instantaneous fluctuations of the voltage V 22 

We finally note that in the experiment of Ref. [13] 
the circuit corrections have not been observed, which 
we explain by a very low impedance of the circuit, 
Gs/Gtot <C 1, operating in the weak damping regime 
with Q = Gcjpi/Gtot = 2.5. Indeed, in this case 
f = {Gs/2G tot )eV s (assuming T = 0), therefore 
the circuit corrections are suppressed by a small fac- 
tor Gs/Gtot- On the other hand, the first term in Eq. 
(11a) is So = AU/T and the second term can be esti- 
mated as (Gtot/Gs)(7lwd/Q)So, where TZ wd = fruj p \/eV s 
in this case. Therefore the relative contribution of the 
third cumulant to the total action increases by the fac- 
tor G 'tot I 'Gs, which makes it favorable to use a low- 
impedance circuit. 

In conclusion, we have formally solved Kramers' prob- 
lem of noise-activated escape from a metastable state 
for Markovian non-equilibrium noise beyond the Gaus- 
sian approximation. We have applied this result to the 
Josephson junction threshold detector and evaluated the 
third current cumulant contribution to the escape rate. 
The back action of the electrical circuit on the measured 
noise is entirely a non-equilibrium effect, which leads to 
"cascade" corrections for the third cumulant. It would 
be interesting to apply our theory to the dynamics of the 
Josephson bifurcation amplifier. 23 
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APPENDIX A 

In this Appendix, we give the details of the derivation 
of Eq. (11). Starting with Eq. (8), we expand the gencr- 
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ating function in the counting variable to third order, 

H = (I S - J5)A in 0+((4))( Ain 0) 2 /2+((/|))(A in 0) 3 /3! + . . . 

(Al) 

Taking the assumptions in the text of an ohmic system 
conductor, we replace (Ig — Jg) = —GgV. The noise 
produced by the system conductor depends on the dy- 
namically changing voltage drop Vg — V (see Fig. 1), so 
we also expand the second current cumulant in the volt- 
age V, namely «j£» = «j£» - d Vs {{I 2 s ))V + . . ., where 
we recall that Vg is the voltage across the system, and 
the coefficients above are evaluated at V = 0. Using the 
Josephson relation V = ^o'P, and averaging (Al) over 
the oscillation period at fixed energy, we now need to 
find the instanton line and solve the equation 







-G s $o(^A in + ((/!)) (4> 2 ) E \?J2 (A2) 
^d Vs (I 2 s))^)E\t/2 + ((/f))^ 3 ) £ A 3 n /3! + . . . 



for X in (E). 

As argued in the text, the instanton solution X- ln (E) is 
small for the J J, so the series may be inverted to find to 
third order in fkj p \/T e g, 



A, 



T eff QT^Gg(^) E ' 



(A3) 



Inserting this instanton solution into the action, we ob- 
tain for the activation rate 

iogr = -^ + ^»- fdE^. (A4) 

However, to leading order in damping, the average 
(0 3 )_e vanishes because 4> is an °dd function of time 



over the interval of one period. Therefore, we must 
evaluate this average to next order in the variation of 
energy over one period. Recalling that the action is 
S = J dt'[-\ E E/<& + H{\e4>)], the time dependence 
of the energy E may be found by taking the functional 
derivative 5S/5\e to find E = $od\ E H. Inserting the 
instanton solution Ai n = &o/T c s, we obtain 



E = $^G S 
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(A5) 



so the energy is increasing. 

At a given value of <f>, the linear correction to 4> 2 due 
to energy dissipation is given by 5<fi 2 — (2/m)SE. There- 
fore, we use dt<j> — d<fi to write 



Jo 



mT 



d(/)5E, (A6) 



where the integral starts from the turning point nearest 
to the top of the barrier, <p = 4>o- This is done because 
the escape of the JJ into the running state happens at 
the top of the barrier, so the last part of the trajectory 
before escape must be fully accounted for. 

Integrating by parts, § dcf>5E = (fiSE^ — §dt4>E, 
we keep only the change in energy, and neglect the 
change in the turning point to this order. Using the re- 
sult (A5), these approximations imply that § d<j)5E = 
$§Gs § dt(4>o — <t>)<i> 2 - Collecting results, we find 



2*§G S 



((^o - 4>)4?) 



m 



(A7) 



Recalling the definition m = &qC, and inserting (A7) 
into (A4), we recover our main result (11). 
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